™hD 


HD-A134  831 


UNCLASSIFIED 


POLARIZATION  PATTERNS  OF  PI  2  MAGNETIC  PULSATIONS  AND 
THE  SUBSTORM  CURRENT  WEDGECU)  BOSTON  UNIV  MA  DEPT  OF 
ASTRONOMV  M  LESTER  ET  AL.  01  OCT  83  AFGL-TR-83-0302 
F19628-81-K-0003  F/G  8/14 


1/1 


BUC  FILE  COPY 


JOURNAL  OF  GEOPHYSICAL  RESEARCH.  VOL.  88,  NO.  AIO.  PAGES  7958-7966.  OCTOBER  I.  1983 


M 


AFG  T, -  T"  -  '—  0302 

Polarization  Patterns  of  Pi  2  Magnetic  Pulsations  and  the  Substorm  Current 

Wedge 

Mark  Lester  and  W.  Jeffrey  Hughes 

Department  of  Astronomy .  Boston  University 

Howard  J.  Singer 


m 


r  *  •  <  > 

^bLECTE 


Air  Force  Geophysics  Laboratory 

The  results  of  an  analysis  of  16  Pi  2  pulsations  observed  with  an  extended  east-west  chain  of  mid¬ 
latitude.  ground-based  magnetometers  are  reported.  The  events  were  chosen  such  that  the  center  of 
the  substorm  current  system,  defined  by  using  the  midlatitude  bays  associated  with  these  pulsations, 
was  within  the  longitudinal  extent  of  the  station  chain.  A  current  wedge  model  is  also  used  to  interpret 
the  observed  polarization  pattern  of  the  Pi  2  pulsations.  This  pattern  is  used  to  locate  the  center  of  the 
Pi  2  current  system.  The  centers  of  the  Pi  2  and  substorm  current  systems  occur  at  the  same  meridian 
for  only  -65%  of  the  events.  This  result  suggests  that  the  Pi  2  current  system  and  the  substorm  current 
wedge  are  not  always  the  same.  The  longitudinal  extent  of  the  Pi  2  current  system  estimated  from  the 
Pi  2  polarization  pattern  averages  ~90°  or  -6  hours  local  time  for  the  events  in  the  study.  The  sense  of 
ellipticity  of  the  waves,  anticlockwise  looking  down  the  held  line  in  the  northern  hemisphere,  agrees 
with  previous  results  from  these  latitudes.  Estimates  of  the  wave-phase  difference  between  stations 
show  that,  in  general,  the  eastern  station  of  a  station  pair  leads  the  western  one  at  all  local  times  for 
both  H  and  D  components.  When  plotted  in  a  substorm  coordinate  system  based  on  the  mid-latitude 
bay.  the  phase  difference  per  degree  of  longitude  shows  a  tendency  to  decrease  in  the  eastern  portion 
of  the  current  wedge.  This  longitudinal  pattern  of  phase  difference  is  consistent  with  the  eastern, 
downward.  Held-aligned  currents  being  less  localized  than  the  western,  upward.  Held-aligned  currents. 
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I .  Introduction 

The  appearance  of  Pi  2  magnetic  pulsations  at  the  onset  of 
substorms  has  been  recognized  for  many  years  [Saito,  1961  ]. 
Recently,  a  consensus  [Rostoker  et  at.,  1980]  included  at 
least  one  Pi  2  pulsation  event  as  an  integral  part  of  a 
substorm,  and  the  Pi  2  signature  has  been  used  to  time 
substorm  onsets  [ Rostoker  and  Olson,  1979,  and  references 
therein].  Because  of  its  association  with  substorm  onset,  the 
study  of  this  class  of  pulsation  is  important  for  the  under¬ 
standing  of  substorm  dynamics.  In  this  paper  we  report 
observations  of  Pi  2  pulsations  made  with  the  AFGL  (Air 
Force  Geophysics  Laboratory)  east-west  chain  of  mid-lati¬ 
tude  magnetometers  [ Knecht  et  a!.,  1979]  to  investigate  the 
relationship  between  the  waves  and  the  substorm  current 
system  [McPherron  et  at.,  1973;  Baumjohann  et  al.,  1981]. 

The  AFGL  magnetometer  chain  is  particularly  well  suited 
for  the  study  of  Pi  2  pulsations  for  several  reasons.  The  chain 
extends  over  4  hours  magnetic  local  time.  This  makes  it 
useful  for  studying  the  three-dimensional  current  systems 
associated  with  Pi  2  pulsations.  Also,  since  Pi  2  pulsations 
can  often  be  observed  over  a  large  portion  of  the  nightside  on 
the  ground  [ Saito  et  al.,  1976;  Sutcliffe,  1980;  Singer  et  al., 
1983]  the  chain  can  be  used  to  study  the  longitudinal 
variations  in  the  wave  structure.  Finally,  the  mid-latitude 
location  of  the  stations  is  a  region  where  the  Pi  2  signal  is 
relatively  free  from  contamination  caused  by  currents  in  the 
auroral  ionosphere  [Samson,  1982],  and  consequently  the 
mid-latitude  signal  provides  clearer  results  than  the  auroral 
zone  signature  [Samson  and  Harrold,  1983], 

The  three-dimensional  substorm  current  wedge  has  been 
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described  by  McPherron  et  al.  [1973],  who  locate  the  center 
of  the  wedge  at  the  longitude  where  there  is  a  change  in  sign 
of  the  east-west  magnetic  bay.  If  it  is  assumed  that  the  Pi  2 
pulsation  results  from  an  oscillation  of  this  current  system 
[Bostrom,  1972],  then  the  mid-latitude  Pi  2  horizontal  polar¬ 
ization  pattern  should  exhibit  a  predictable  longitudinal 
pattern.  The  predicted  polarization  pattern  provides  another 
technique  for  locating  the  center  of  the  wedge  [ Knecht , 
1981].  Comparing  the  two  techniques,  we  find  that  the 
methods  agree  in  only  approximately  two  thirds  of  the  cases. 
Little  previous  work  has  made  use  of  the  substorm  current 
wedge  to  order  the  Pi  2  polarization  [Knecht,  1981;  Samson 
and  Harrold,  1983],  although  in  the  past,  several  apparent 
local  time  variations  have  been  reported  [Bjornsson  et  al., 
1971;  Fukunishi,  1975;  Baranskiy  et  al.,  1980;  Stuart  and 
Baranskiy,  1982], 

The  phase  differences  between  stations  for  Pi  2  pulsations 
are  examined.  Although  there  have  been  some  differences  in 
the  results  of  previous  studies  [ Herron ,  1966;  Mier-Jedr- 
zejowicz  and  Southwood,  1979;  Baranskiy  et  al.,  1980],  the 
predominant  observation  is  that  eastern  stations  lead  in 
phase.  This  result  has  been  interpreted  as  western  phase 
propagation  on  the  entire  nightside.  Since  no  physical  inter¬ 
pretation  consistent  with  substorm  models  has  been  devel¬ 
oped,  we  concluded  that  a  further  study  of  phase  variations 
was  warranted.  Our  results  confirm  that  eastern  stations  lead 
in  phase,  but  that  the  phase  variation  is  not  uniform  across 
the  substorm  current  wedge. 

The  following  section  describes  the  magnetometer  array 
and  data  set  used.  We  then  describe  the  model  of  the 
substorm  current  wedge  and  the  associated,  predicted  Pi  2 
polarization  pattern.  Examples  of  Pi  2  pulsations,  their 
polarization  pattern,  and  the  related  bay  structure  are  then 
given.  The  polarization  pattern  is  described  for  16  events 
where  it  is  possible  to  locate  the  center  of  the  substorm 
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Fig.  I.  Map  showing  the  locations  of  the  AFGL  network  sta¬ 
tions.  The  five  northern  stations,  extending  over  4  hours  local  time, 
are  at  —55°  corrected  geomagnetic  latitude:  the  two  southern 
stations,  which  are  not  used  in  this  study,  are  at  -40°. 

current  system  between  the  two  extreme  meridians  of  the 
network.  Finally  our  results  are  discussed  in  the  context  of 
previous  studies  and  the  three-dimensional  current  system 
model. 

2.  The  Data  Set 

The  source  of  data  for  the  study  is  the  AFGL  magnetome¬ 
ter  network  [Knecht  et  al.,  1979].  The  geographic  and 
corrected  geomagnetic  coordinates  of  the  seven  AFGL 
stations  are  given  in  Table  1 .  Five  of  these  stations  are  at 
—55°  corrected  geomagnetic  latitude,  and  the  remaining  two 
stations  are  at  —40°  (Figure  1).  The  latitudes  of  the  stations 
place  them  all  on  field  lines  which  thread  the  plasmasphere, 
although  it  is  possible  that  during  extremely  disturbed  mag¬ 
netic  conditions  the  plasmapause  may  move  far  enough 
radially  inward  to  place  the  field  lines  of  the  northern 
stations  in  the  plasma  trough.  Data  from  the  two  southern 
stations  are  not  used  in  this  study.  The  total  longitudinal 


separation  of  the  five  northern  stations  is  62.3°  corrected 
geomagnetic  longitude,  which  corresponds  to  slightly  more 
than  4  hours  magnetic  local  time.  The  longitudinal  separa¬ 
tions  between  adjacent  stations  at  the  same  latitude  are  given 
in  Table  1. 

At  each  station  a  three-component  flux  gate  magnetometer 
samples  the  field  once  a  second  with  an  amplitude  resolution 
of  -0.06  nT.  The  three  components  are  H(  +  ve  geomagnetic 
north),  D  (+ve  geomagnetic  east),  and  Z  (  +  ve  vertically 
down).  This  study  uses  5-s  averages.  The  time  interval  from 
which  the  data  were  taken  is  March  14-19,  1978.  During  this 
interval,  we  identified  more  than  40  Pi  2  pulsations  at  one  or 
more  stations. 

3.  The  Substorm  Current  Wedge 
Coordinate  System 

In  many  previous  investigations  of  Pi  2  pulsations,  local 
time  has  been  used  to  order  Pi  2  parameters  such  as  the 
horizontal  wave  ellipticity  [e.g.,  Fukunishi ,  1975;  Lester  and 
On,  1981],  the  orientation  of  the  major  axis  of  the  horizontal 
ellipse  [e.g.,  Baranskiy  et  al.,  1980:  Stuart  and  Baranskiy, 
1982],  and  the  phase  difference  between  stations  [e.g.. 
Herron,  1966;  Mier-Jedrzejowicz  and  Southwood,  1979], 
with  occasional  differences  between  studies  being  reported. 
Since  Pi  2  pulsations  are  associated  with  substorms,  it  seems 
natural  to  use  the  position  of  the  stations  relative  to  some 
substorm  feature  to  better  organize  the  data.  Features  such 
as  the  auroral  breakup  region  [Kuwashima  and  Saito,  1981] 
or  the  auroral  electrojet  [Rostoker  and  Samson,  1981;  Sam¬ 
son  and  Rostoker,  1983]  would  be  useful,  but  no  method  of 
reliably  locating  them  with  mid-latitude  data  is  available. 
However,  with  an  east-west  chain  of  mid-latitude  magne¬ 
tometers,  one  can  locate  the  position  of  the  substorm  current 
wedge  by  using  a  model  described  by  McPherron  et  al. 
[1973].  The  local  time  of  isolated  substorm  onsets  can  vary 
substantially,  and  there  is  usually  a  westward  motion  during 
a  sequence  of  substorm  intensifications  [  Wiens  and  Ros¬ 
toker,  1975;  Knecht,  1981].  Consequently,  organizing  the 
Pi  2  parameters  by  using  the  location  of  the  wedge  relative  to 
a  station  has  clear  advantages  over  the  use  of  local  time. 

The  simple  current  system  discussed  by  McPherron  et  al. 
[1973]  and  developed  further  by  Clatter  and  McPherron 
[1974a]  is  used  as  the  model  and  is  shown  in  Figure  2.  At  the 
onset  of  the  expansion  phase  of  the  substorm  the  crosstail 
current  is  diverted  down  the  magnetic  field  lines.  Current 
then  flows  in  the  ionosphere  as  the  westward  electrojet  and 
returns  to  the  tail  along  the  magnetic  field  lines.  The  pertur¬ 
bation  tail  current  can  be  represented  by  an  equivalent 
eastward  current,  which  completes  the  three-dimensional 
current  wedge.  Although  this  current  system  is  an  oversim¬ 
plification  several  authors  [e.g..  McPherron  et  al.,  1973; 


TABLE  1.  Name  and  Coordinates  of  AFGL  Magnetometer  Network  Stations 


Geographic  Corrected  Geomagnetic 

East  North  East  North  Station 

Station  Name  Longitude  Latitude  Longitude  Latitude  Separation 

NEW,  Newport.  WA  242.9° 

RPC.  Rapid  City,  SD 
CDS,  Camp  Douglas,  WI 
MCL,  Mt.  Clemens.  Ml 
SUB.  Sudbury,  MA 
LOC,  Lompoc,  CA 
TPA,  Tampa,  FL 


48.3°  299.6°  St 

44.2  3 
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Fig.  2.  The  upper  panel  shows  schematically  the  model  of  the 
substorm  current  wedge  which  is  used  to  position  the  stations 
relative  to  the  center  of  the  substorm  wedge.  The  bottom  panel 
(right)  shows  an  equivalent  current  system  and  model  parameters. 
The  bottom  panel  (left)  summarizes  the  calculated  sense  and  magni¬ 
tude  of  the  mid-latitude  magnetic  bay  expected  after  substorm  onset 
[from  Clauer  and  McPherron,  1974a], 


Clauer  and  MiPherron,  1974a]  have  shown  that  it  accurately 
describes  mid-latitude  bay  variations.  The  mid-latitude  sig¬ 
nature  of  this  wedge  is  a  positive  perturbation  in  the  north- 
south  magnetic  component  symmetric  about  the  central 
meridian  of  the  current  system  and  a  perturbation  in  the  east- 
west  magnetic  component  that  is  positive  to  the  west  and 
negative  to  the  east  of  the  central  meridian. 

Figure  3  shows  schematically  the  field-aligned  and  iono¬ 
spheric  legs  of  the  substorm  current  wedge  and  the  magnetic 
variations  caused  by  the  two  field-aligned  currents.  The 


field-aligned  currents  will,  ideally,  only  affect  the  horizontal 
magnetic  components  at  mid-latitudes.  If  the  Pi  2  pulsation  is 
directly  associated  with  the  substorm  current  system,  then  it 
will  have  the  predicted  polarization  pattern  across  the  cur¬ 
rent  wedge  shown  in  Figure  3.  This  pattern  can  be  tested  by 
using  events  where  the  center  of  the  current  wedge,  located 
by  the  bay  current  system,  is  between  the  two  longitudinal 
extremes  of  the  magnetometer  chain. 

In  order  to  quantify  the  bay  structure  for  each  Pi  2  event 
we  measured  both  the  H  and  D  values  at  the  five  northern 
stations  at  the  start  of  the  Pi  2  pulsation  and  at  a  time  when 
no  more  wave  activity  was  observed.  We  defined  the  differ¬ 
ence  between  the  two  values  as  A H  and  AD.  The  time 
difference  was  typically  between  10  and  15  minutes.  The 
central  meridian  is  where  ID  =  0.  The  H  perturbation 
provides  a  consistency  check,  i.e..  IH  should  be  positive 
and  maximize  near  the  location  of  ID  =  0.  Of  the  original  40 
events.  32  had  such  a  bay  structure,  and  for  16  of  these  32 
the  location  of  ID  =  0  was  between  the  NEW  and  SUB 
meridians. 

The  longitude  of  ID  =  0  was  calculated  by  assuming  a 
linear  change  in  ID  between  the  two  stations  on  either  side 
of  ID  =  0.  Each  station  was  located  in  longitude,  relative  to 
the  meridian  where  AD  =  0.  by  calculating  the  longitude 
difference:  Along.  Along  was  positive  (negative)  if  the  station 
was  to  the  east  (west)  of  the  AD  =  0  longitude. 

4.  Examples  of  P)  2  Pulsations  and  the 
Substorm  Current  System 

Before  discussing  the  results  of  all  16  events  statistically, 
two  events  are  used  to  illustrate  the  analysis  procedure.  Pi  2 
polarization  parameters  were  determined  both  by  plotting 
hodograms  and  by  using  spectral  analysis  techniques.  Power 
spectra  were  computed  for  each  field  component  at  each 
station.  The  peak  frequency  in  the  range  5-25  mHz  (200-40 
s)  was  measured  in  the  H  and  D  components  at  each  station. 
The  dominant  frequency  at  a  majority  of  stations  in  each 
magnetic  component  was  chosen  as  the  overall  dominant 
frequency  Polarization  parameters  of  the  Pi  2  pulsation, 
percent  polarization,  ellipticity,  and  azimuth  were  measured 
at  the  dominant  frequency  by  using  standard  spectral  analy¬ 
sis  methods  [see,  e.g.,  Arthur  et  al.,  1976],  The  azimuth,  or 
orientation  of  the  major  axis,  is  positive  (negative)  in  the 
northeast  (northwest)  quadrant.  The  ellipticity  is  0  when 
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Fig.  4.  A  Pi  2  pulsation  observed  by  the  AFGL  network  in  the  interval  0247-0302  UT  on  March  14.  1978.  In  the 
upper  panels  the  unfiltered  H  and  P  components  show  the  Pi  2  pulsation  and  the  magnetic  bay.  The  lower  panels  show 
the  hodograms  measured  during  the  interval  0253-0257  UT  for  each  station.  The  dashed  lines  in  some  of  the  hodograms 
represent  data  gaps.  The  arrow  at  the  side  of  each  hodogram  represents  the  sense  of  ellipticity.  The  hodograms  at  the 
four  eastern  stations  compare  favorably  with  the  model  Pi  2  azimuth  pattern  inside  the  wedge  illustrated  in  Figure  3.  The 
polarization  at  station  NEW  is  that  expected  for  a  station  to  the  west  of  the  wedge. 


linear  and  I  when  circular.  A  positive  (negative)  ellipticity 
indicates  a  clockwise  (anticlockwise)  sense  when  looking 
down  the  magnetic  field  line  in  the  northern  hemisphere. 
Temporal  variations  of  ellipticity  and  azimuth  were  exam¬ 
ined  by  plotting  hodograms.  The  two  methods,  hodograms 
and  spectral  wave  analysis,  produced  similar  results  for  the 
ellipticity  and  azimuth  in  all  cases. 

Figures  4  and  5  show  unfiltered  data  for  the  time  intervals 
0247-0302  UT  and  0555-0610  UT  on  March  14.  1978.  each  of 
which  contains  a  Pi  2  pulsation.  The  lower  panels  show  the 
hodograms  for  each  event.  During  the  earlier  event  (Figure 


4).  the  hodogram  at  NEW,  the  station  with  the  smallest  wave 
amplitude,  does  not  fit  the  predicted  pattern  inside  the 
current  wedge  (Figure  3).  However,  it  is  consistent  with  the 
pattern  predicted  to  the  west  of  the  western  field-aligned 
current  part  of  the  wedge  (not  shown  in  Figure  3).  At  the 
other  stations  the  polarization  pattern  is  similar  to  the 
predicted  one  and  points  to  a  location  between  MCL  and 
SUB.  There  is  also  some  time  variation  of  the  azimuth  at 
SUB.  but  the  average  azimuth  fits  the  predicted  pattern.  For 
the  later  event,  shown  in  Figure  5,  the  predicted  pattern  of 
the  azimuth  (Figure  3)  is  reproduced  almost  exactly. 
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Fig  5.  Same  as  Figure  4  for  a  Pi  2  pulsation  in  the  interval  O'SS-OMO  I  I  on  March  14.  I97g  When  no  arrow  is 
present  the  elliplicily  is  essentially  linear  The  holograms  between  0X59-116(1'  IT  from  all  five  stations  compare 
favorably  with  the  model  shown  in  Figure  3 
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Fig.  6.  The  calculated  values  of  A D  (upper  panels)  and  A H 
(middle  panels)  for  the  Pi  2  events  in  Figures  4  and  5  are  shown.  The 
longitudinal  variation  of  AD  and  A H  can  be  compared  with  the 
predicted  mid-latitude  signature  in  Figure  2.  The  location  where  AD 
=  0.  corresponding  to  the  center  of  the  wedge,  is  marked  by  a 
dashed  line  (C).  The  position  of  line  C  can  be  compared  with  the 
location  where  the  major  axis  of  the  polarization  ellipses  (lower 
panels)  points  for  the  two  events. 


The  bay  structures  (Figures  4  and  5,  upper  panels)  show 
the  schematic  variation  of  the  bays  discussed  in  the  previous 
section.  The  change  in  sign  of  the  D  bay  occurs  between 
NEW  and  RPC  for  the  0247  UT  event  and  between  RPC  and 
CDS  for  the  0555  UT  event.  Note  also  that,  except  for  MCL 
and  SUB,  for  the  first  event  the  H  bays  are  positive.  The  D 
bay  reaches  the  most  negative  value  for  MCL,  which  sug¬ 
gests  that  MCL  and  SUB  are  located  at  longitudes  outside  of 
the  wedge. 

The  bay  structure  variation  is  further  illustrated  in  Figure 
6  by  using  the  values  of  A//  and  A D.  For  both  events  the 
upper  panel  shows  the  A D  values  for  each  station,  and  the 
middle  panel  shows  the  A H  values.  The  position  where 
A D  =  0  is  estimated,  and  it  is  marked  by  a  dashed  line 
indicating  the  wedge  center  (C).  For  the  event  at  0247  UT  the 
line  C  passes  between  NEW  and  RPC,  but  the  polarization 
pattern  (bottom  left  panel),  although  similar  to  the  predicted 
one.  points  between  MCL  and  SUB,  as  mentioned  earlier. 
This  type  of  event,  where  the  two  methods  of  identifying  the 
center  of  the  wedge  disagree,  was  found  on  five  occasions. 
The  size  of  the  disagreement  given  in  the  example  is  the  most 
extreme,  but  in  all  five  cases  the  difference  is  more  than  a 
simple  displacement  to  the  wrong  side  of  a  station.  The 
second  event  at  0555  UT  shows  that  the  center  of  the  wedge 
identified  by  AD  =  0  agrees  with  the  center  found  by  the 
polarization  pattern.  This  type  of  event  was  found  on  1 1 
occasions. 

5.  Polarization  Variations  Relative  to  the 
Substorm  Current  Wedge 

Having  shown  two  examples  of  the  polarization  patterns 
of  Pi  2  pulsations,  we  now  present  the  polarization  parame¬ 
ters  for  the  16  events  where  the  center  of  the  current  wedge 
was  located  between  the  NEW  and  SUB  meridians.  The 
azimuth  and  ellipticity  were  computed  by  using  the  spectral 


analysis  method  with  2.67  degrees  of  freedom.  The  conven¬ 
tions  for  azimuth  and  ellipticity  were  described  in  the 
previous  section.  Only  values  for  which  the  percent  polariza¬ 
tion  exceeded  959f  were  considered.  If  the  modulus  of  the 
ellipticity  |f|  >  0.80.  then  the  azimuth  value  was  not  consid¬ 
ered,  since  the  pulsation  is  nearly  circularly  polarized  in 
such  cases,  and  determining  azimuth  is  imprecise  |cf.  Green 
and  Hamilton.  1978,  1981). 

The  variation  of  azimuth  for  the  16  events  where  the 
center  of  the  wedge  is  identified  between  NEW  and  SUB  is 
shown  in  Figure  7.  The  azimuth  at  each  station  for  each 
event  is  plotted  in  the  substorm  current  wedge  coordinate 
system  where  Along  =  0  is  at  the  location  of  AD  =  0.  Only 
values  that  satisfied  the  selection  criteria  described  in  the 
above  paragraph  are  plotted.  The  events  where  the  two 
methods  for  defining  the  center  of  the  substorm  current 
wedge  agree  (dots)  show  the  expected  change  of  polarization 
across  the  wedge.  For  purposes  of  comparison  the  dashed 
line  represents  the  Pi  2  azimuths  expected  from  a  substorm 
current  wedge  with  a  longitudinal  extent  of  90°  symmetric 
about  the  central  meridian.  The  polarization  azimuths  for 
those  events  where  the  two  methods  disagree  are  indicated 
by  crosses.  The  crosses  corresponding  to  each  event  have 
been  joined  by  straight  lines.  The  remarkable  feature  is  that 
these  events  also  show  a  pattern  similar  to  the  dots  but 
displaced  relative  to  the  wedge  center.  The  displacement  is 
not  consistently  to  one  side  or  the  other  of  the  dashed  line. 

Azimuth 

West 


North 


East 

A  Long 

Fig.  7.  Azimuth  versus  Along  for  (he  16  Pi  2  events  discussed  in 
the  study.  Tick  marks  on  (he  azimuth  axis  are  every  30”.  The  dots 
represent  the  azimuth  of  the  events  where  the  two  methods  for 
locating  the  center  of  the  current  system  agree.  The  crosses  repre¬ 
sent  the  azimuth  for  the  events  where  the  two  methods  disagree. 
The  dashed  line  represents  the  variation  expected  for  a  model 
current  system  6  hours  wide  and  symmetric  about  the  central 
meridian.  The  error  in  the  measured  azimuth  value  is  HI  . 
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Fig.  8.  Ellipticity  versus  Along  for  the  16  Pi  2  events  discussed 
in  the  study.  Positive  (negative)  ellipticity  indicates  a  clockwise 
(anticlockwise)  sense  of  rotation  looking  down  the  field  line  in  the 
northern  hemisphere.  The  center  of  the  wedge  is  located  at  Along  -  0. 

indicating  the  center  of  the  Pi  2  current  system  can  be  either 
to  the  east  or  west  of  the  substorm  current  wedge.  It  should 
be  noted  that,  although  we  do  not  present  the  results,  we 
found  no  systematic  variation  of  azimuth  with  local  time. 

In  the  substorm  current  wedge  model  (Figure  3)  the 
position  where  the  azimuth  of  the  magnetic  perturbation  is 
directed  east-west  locates  the  longitudes  of  the  field-aligned 
current  parts  of  the  wedge.  The  difference  in  longitude 
between  the  field-aligned  current  portions  can  be  regarded  as 
the  effective  width  of  the  wedge.  Since  it  is  unusual  for  the 
entire  width  of  the  wedge  to  be  located  within  the  two 
longitudinal  extremes  of  the  network,  the  wedge  width  was 
estimated  from  the  slope  of  the  lines  connecting  the  dots  or 
crosses  for  each  individual  event.  This  slope  varies  from  - 1° 
to  -2.5°  azimuth/degree  longitude,  with  an  average  of  ~2° 
azimuth/degree  longitude.  This  implies  a  variation  in  the 
wedge  width  of  between  -4  and  -12  hours  and  an  average  of 
—6  hours.  These  results  are  in  reasonable  agreement  with 
the  observations  of  Clatter  and  McPherron  (I974/>|.  who 
found  the  extent  of  the  current  system  to  vary  between  2.5 
and  10  hours,  with  an  average  of  just  over  4  hours. 

The  ellipticity  variation  with  respect  to  Along  is  shown  in 
Figure  8,  where  all  but  two  values  are  negative.  The  predom¬ 
inance  of  anticlockwise  values  agrees  with  other  results  for 
stations  at  these  latitudes  [e.g..  Fitkinishi.  1975;  Mier-Jedr¬ 
zejowicz  and  Southwood.  1979;  Lester  and  Orr.  1 98 1 1 .  Apart 
from  this,  no  ordering  appears  «n  the  sense  of  polarization 
with  respect  to  Along.  Also,  we  found  no  systematic  local 
time  variation  of  ellipticity. 

6.  Phase  Differences  of  Pi  2  Pulsations  Between 
Stations 

The  longitudinal  variation  of  pulsation  phase  has  been 
discussed  in  the  past  in  terms  of  the  rn  number  |c.g..  Green. 


1976:  Mier-Jedrzejowicz.  and  Southwood.  I979|.  The  in  num¬ 
ber  is  obtained  by  dividing  the  phase  difference  between  a 
pair  of  stations  by  the  station  longitudinal  separation.  It  has 
then  been  interpreted  in  terms  of  an  apparent  azimuthal 
wave  number.  However,  since  the  azimuthal  propagation  of 
the  Pi  2  signal  is  not  clear,  we  shall  not  interpret  the  ni 
number  in  this  way  but  simply  consider  it  as  the  phase 
difference  per  degree  of  longitude. 

The  phase  differences  between  all  possible  pairs  of  sta¬ 
tions  were  calculated  by  using  a  cross-spectral  analysis 
technique  | Hitultes  el  at..  I978|.  Only  the  phase  differences 
between  adjacent  station  pairs  are  used  here.  The  9 5r4 
confidence  limits  in  the  phase  difference  were  calculated  by 
using  the  coherency  value  and  the  number  of  degrees  of 
freedom  associated  with  the  smoothing  parameter  used  in 
the  spectral  computation  | Green.  I976|.  We  calculated  the  m 
number  from  the  calculated  phase  difference  and  an  error  in 
m.  Am.  from  the  confidence  limits  in  the  phase  difference. 

Figure  9  shows  the  m  values  for  the  1 1  events  where  the 
two  methods  of  identifying  the  center  of  the  substorm 
current  wedge  agree.  The  error  bar.  representing  Am.  is 
plotted  if  Am  >  |m|.  i.e..  there  is  ambiguity  in  the  sign  of  m. 
We  included  only  those  points  for  which  Am  <  1.5.  This 
value  was  chosen  as  a  compromise  between  having  sufficient 
number  of  points  and  sufficient  resolution  in  m.  The  maxi¬ 
mum  error  allowed.  Am  -  1.5.  is  also  shown.  The  // 
component  (Figure  9«)  shows  a  weak  tendency  for  im|  to 
decrease  from  west  to  east  across  the  wedge.  The  D  compo¬ 
nent  (Figure  9b)  shows  a  similar  tendency  for  the  values  to 
the  east  of  the  center  to  decrease  with  increasing  Along.  The 
larger  jm|  values  are  west  of  center  of  the  wedge  for  the  // 
component  and  close  to  the  center  for  the  D  component. 

The  m  values  are  predominantly  negative,  irrespective  of 
the  location  of  the  center  of  the  substorm  current  wedge  and 
local  time  (results  not  presented).  This  agrees  with  the  work 
of  Mier-Jedrzejowicz  and  Southwood  1 1979|.  Warner  1 1979). 
and  Green  and  Stuart  [  19791.  Mier-Jedrzejowicz  and  South- 
wood  did  not  use  individual  events,  but  power  levels  in  the 
Pi  2  frequency  band  averaged  over  I  hour  time  intervals. 
Warner,  however,  studied  individual  Pi  2  events  observed 
with  the  same  network  as  Mier-Jedrzejowicz  and  South- 
wood.  Baranskiy  el  al.  [  1980)  reported  that  there  was  a  slight 
tendency  for  m  to  change  sign  near  midnight.  They  have 
shown  the  H  component  of  the  western  station  leading 
before  midnight  and  the  l)  component  of  the  western  station 
leading  after  midnight.  Neither  reversal  is  seen  in  our 
results.  Also,  there  is  no  sign  of  the  reversal  reported  by 
Herron  1 1966|.  who  used  data  from  a  similar  longitude  sector 
and  latitude  as  the  AFGL  network. 

7.  DtsuussioN 

For  each  of  the  16  events  studied  in  detail  the  azimuths  of 
the  major  axis  of  the  polarization  ellipses  agree  with  the 
pattern  predicted  by  a  model  three-dimensional  current 
system  (Figure  2).  The  pattern  was  better  organized  by  the 
substorm  current  wedge  and  not  by  local  time.  Measure¬ 
ments  made  by  Green  and  Stuart  |  1979;  personal  communi¬ 
cation.  I980|  and  Baranskiy  el  al.  1 198()|  support  our  results. 
Two  examples  of  the  azimuth  of  (he  major  axis  of  the 
polarization  ellipse  |see  Figure  5.  Haranskiv  et  al..  I98(t| 
illustrate  the  pattern  predicted  by  a  current  wed<’e  model. 
Furthermore,  using  many  events.  Baranskiy  et  al.  did  not 
find  a  systematic  variation  with  local  time  (see  their  Figure 
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Fig.  9.  The  phase  difference  per  degree  of  longitude,  m  value, 
between  pairs  of  stations  plotted  against  the  Along  of  the  midpoints 
of  the  two  stations  for  H  (9a)  and  D  (9b);  m  is  calculated  by  using 
only  adjacent  station  pairs  for  the  1 1  events  where  the  two  methods 
of  locating  the  center  of  the  current  system  agree.  Positive  (nega¬ 
tive)  m  indicates  the  western  (eastern)  station  leads  in  phase.  Error 
bars  are  marked  if  A m  >  |m|  and  the  maximum  error  in  any  value 
plotted  is  1.5. 

6c).  However,  there  have  been  some  studies  which  suggest  a 
local  time  organization  of  the  Pi  2  polarization  pattern  [e.g.. 
Bjornsson  et  al.,  1971;  Fukinishi,  1975;  Stuart  and  Barans- 
kiy,  1 982| .  Bjornsson  et  al.  [1971],  using  averages  of  Pi  2 
azimuths  determined  in  I  hour  local  time  bins  at  stations 
between  40°  and  50°  magnetic  latitude,  found  a  pattern 
similar  to  that  in  Figure  3.  but  pointing  to  the  2200  LT 
meridian.  These  correlations  between  Pi  2  azimuths  and 
local  time  would  be  expected  if  the  centers  of  the  current 
systems  associated  with  the  pulsations  used  in  a  study  had  a 
small  local  time  spread.  We  are  unable  to  test  this  explana¬ 
tion  because  the  studies  showing  a  local  time  organization  of 
the  Pi  2  azimuths  do  not  indicate  the  location  of  the  wedge 
centers.  In  this  paper  the  centers  of  the  current  systems  were 
found  to  range  from  1940  to  0100  LT.  a  range  that  makes  it 


possible  to  distinguish  between  a  local  time  or  current  wedge 
organization  of  the  data. 

The  occasional  disagreement  between  the  center  of  the 
wedge  identified  by  the  two  methods  implies  that  if  the  Pi  2 
pulsation  is  associated  with  a  three-dimensional  current 
system  then  that  current  system  is  not  always  colocated  with 
the  substorm  current  system.  A  study  by  Pashin  et  al.  [  I982| 
appears  to  support  this  result,  as  they  also  conclude  that  the 
wave  and  sulstorm  current  systems  are  offset.  However, 
since  they  were  using  data  from  a  small  region  at  high 
latitudes,  the  exact  relationship  between  their  results  and 
ours  requires  further  study.  Alternatively,  the  Pi  2  could  be 
linked  to  the  substorm  current  system  only  at  onset  (J.  C. 
Samson,  personal  communication.  1982).  After  onset  there 
often  is  a  westward  expansion  of  the  current  system,  which 
will  affect  the  bay  at  mid-latitudes.  In  this  case  the  center 
identified  by  A D  =  0  would  be  to  the  west  of  that  identified 
by  the  Pi  2  polarization  pattern.  This  happens  in  four  out  of 
five  of  the  cases  where  the  two  methods  disagreed. 

As  previously  mentioned,  we  believe  the  wedge  current 
system  discussed  in  this  paper  is  an  oversimplification  of  the 
real  current  systems  [ Baamjohann  et  al.,  1981 1.  There  are 
several  situations  that  could  account  for  discrepancies  be¬ 
tween  observations  and  model.  First,  there  are  likely  to  be 
current  systems  located  outside  the  current  wedge.  Second, 
preexisting  currents  in  the  vicinity  of  the  substorm  current 
wedge  can  confuse  the  bay  pattern.  Third,  the  currents  may 
be  distributed  in  both  latitude  and  longitude.  Nevertheless, 
we  have  demonstrated  that  the  simple  wedge  model  often 
orders  the  data  well. 

In  this  paper  we  are  not  attempting  to  completely  describe 
the  currents  that  are  an  integral  part  of  the  Pi  2  wave.  Since 
Pi  2  polarizations  are  generally  elliptical,  the  current  patterns 
must  move  periodically.  What  we  have  found  is  that  the 
largest  magnetic  perturbations,  those  corresponding  to  the 
major  axes  of  the  polarization  ellipses,  do  fit  a  current 
topologically  similar  to  the  substorm  current  system.  The 
magnetic  signature  seen  at  one  ground  station  describes  the 
integrated  effect  of  currents  flowing  in  a  substantial  region  of 
the  ionosphere.  If  the  ionosphere  conductivity  is  uniform, 
the  ground  signature  results  from  Hall  currents  alone 
[Hughes  and  Southwood,  1976].  However,  at  night  the 
conductivity  is  far  from  uniform;  the  auroral  zone  can  have  a 
substantially  higher  conductivity  than  the  surrounding  iono¬ 
sphere.  In  this  case  the  ground  signature  can  only  be 
determined  by  a  full  Biot-Savart  integration  of  the  iono¬ 
spheric  and  field-aligned  currents.  The  inverse  problem  is 
inherently  nonunique.  A  full  description  of  the  currents 
associated  with  a  Pi  2  pulsation  must  await  further  study. 

The  predominance  of  negative  m  values  in  both  H  and  D 
components  found  in  section  6  agrees  with  Mier-Je- 
drzejowicz  and  Southwood  1 1979]  and  Warner  1 1979|  but  not 
Herron  [1966],  Previously,  nt  values  have  been  interpreted 
as  plane-wave  wave  numbers,  which  in  this  case  implies 
westward  phase  propagation  across  the  entire  nightside.  We 
consider  it  dangerous  to  interpret  m  values  in  this  way.  To 
do  so  implies  a  Pi  2  wave  source  on  the  morningside  of  the 
magnetosphere,  and  we  have  seen  no  substorm  model  that 
can  incorporate  this  interpretation,  since  substorms  are  a 
nightside  phenomenon  with  a  peak  occurrence  between  22 
LT  and  local  midnight. 

Finally,  the  weak  tendency  of  |/»|  to  decrease  to  the  east  of 
the  center  of  the  wedge  is  consistent  with  the  eastern  field- 
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aligned  currents  being  located  over  a  wider  longitudinal 
range  than  the  more  intense  upward  field-aligned  currents  at 
the  western  part  of  the  current  system.  This  result  is 
suggested  by  the  tendency  for  localized  structures  to  be 
associated  with  higher  in  numbers.  Such  a  current  system 
has  been  implied  from  ground-based  magnetometer  mea¬ 
surements  and  STARE  electric  field  measurements  in  Scan¬ 
dinavia  [ Baumjohann  el  of.  1981).  Westward  motion  of  such 
a  current  system  relative  to  the  ground  observations  would 
mean  that  the  phase  change  would  be  larger  at  the  western 
part  of  the  array,  and  the  eastern  stations  would  lead  in 
phase.  Such  motion  of  the  substorm  current  system  does 
seem  to  occur  [e.g..  Wiens  and  Rostoker.  1975). 

8.  Summary  and  Conclusions 

Results  from  a  study  of  Pi  2  pulsations  observed  on  an 
east-wed  chain  of  magnetometers  have  been  presented.  We 
have  concentrated  on  the  relationship  between  magnetic 
bays  and  Pi  2  polarization  characteristics.  We  also  calculated 
the  phase  differences  between  adjacent  station  pairs.  We  can 
summarize  the  results  as  follows: 

1.  The  magnetic  bays  observed  at  mid-latitude  ground 
stations  can  be  used  to  locate  the  center  of  the  substorm 
current  wedge. 

2.  A  predicted  Pi  2  polarization  pattern  based  on  a 
current  wedge  model  is  observed  in  each  of  16  events. 

3.  Two  methods  for  locating  the  center  of  the  substorm 
current  wedge  agreed  in  1 1  out  of  16  cases. 

4.  Using  the  Pi  2  azimuths,  the  effective  width  of  the  Pi  2 
current  system  was  estimated  to  have  an  average  value  of  6 
hours  for  the  16  events. 

5.  There  is  no  systematic  variation  in  ellipticity  with 
either  local  time  or  the  substorm  current  system. 

6.  Wave  phase  differences  between  stations,  or  in,  are 
predominantly  negative  in  both  H  and  D  components. 

7.  Wave  phase  difference  between  stations  tends  to 
decrease  to  the  east  of  the  center  of  the  substorm  current 
wedge,  consistent  with  western  upward  field-aligned  cur¬ 
rents  being  more  localized  than  the  eastern  downward 
currents. 

The  simple  three-dimensional  substorm  current  wedge 
model  orders  Pi  2  polarization  data  well.  The  relationship 
between  the  Pi  2  polarization  and  the  substorm  bay  shows 
that  the  two  current  systems  are  not  always  colocated.  This 
is  important  both  in  terms  of  the  generation  mechanism  of 
the  Pi  2  pulsation  and  also  in  gross  substorm  dynamics. 
Further  study  of  this  property  will  be  carried  out.  and  it  is 
hoped  that  the  results  will  provide  more  information  on  the 
processes  involved  in  the  generation  of  Pi  2  pulsations  and 
their  relationship  with  substorm  onset. 
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